Effect of halo modeling on WIMP exclusion limits 
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WIMP direct detection experiments are just reaching the sensitivity required to detect galactic 
dark matter in the form of neutralinos. Data from these experiments are usually analyzed under the 
simplifying assumption that the Milky Way halo is an isothermal sphere with maxwellian velocity 
distribution. Observations and numerical simulations indicate that galaxy halos are in fact triaxial 
and anisotropic. Furthermore, in the cold dark matter paradigm galactic halos form via the merger 
of smaller subhalos, and at least some residual substructure survives. We examine the effect of 
halo modeling on WIMP exclusion limits, taking into account the detector response. Triaxial and 
anisotropic halo models, with parameters motivated by observations and numerical simulations, 
lead to significant changes which are different for different experiments, while if the local WIMP 
distribution is dominated by small scale clumps then the exclusion limits are changed dramatically. 

98.70.V, 98.80.C 



I. INTRODUCTION 

Arguably the best motivated non-baryonic dark matter 
candidate is the neutralino (the lightest supersymmet- 
ric particle), and current direct detection experiments 
are just reaching the sensitivity required to probe the 
relevant region of parameter space The most strin- 
gent exclusion limits on Weakly Interacting Massive Par- 
ticles (_WIMPs) in general currently come from the Edel- 
weiss [H and Cryogenic Dark Matter Search (CDMS) ex- 
periments with competitive constraints also having 
been produced by Heidelberg-Moscow |^ and IGEX 
The sensitivity to WIMPs will be improved significantly 
in the short term future by the continued operation of 
Edelweiss, and CDMS moving in a low background envi- 
ronment at the Soudan mine and in the longer term 
by, for instance, the planned GENIUS project 

The direct detection event rate, and its energy distri- 
bution, depend crucially on the WIMP speed distribu- 
tion. Data analyzes nearly always assume a standard 
smooth halo model with isotropic maxwellian velocity 
distribution. The change in the expected signal has been 
calculated for various non-standard halo models of vary- 
ing degrees of sophistication |^-pd|. For models which 
are effectively close to maxwellian, while there may be a 
significant change in the annual modulation and angular 
dependence of the signal, the change in the mean (aver- 
aged over time and recoil direction) differential event rate 
is typically small ||]. Models with triaxiality or veloc- 
ity anisotropy may however produce a significant change 
even in the mean differential event rate P,p^. Further- 
more all of the non-standard halo models which have pre- 
viously been considered are essentially smoothFI. N-body 



simulations, however, produce dark matter halos which 
as well as being triaxial with anisotropic velocity distri- 
butions ||l3| , [l^ also contain substructure [|l5|. A number 
of groups have recently investigated the local dark matter 
distribution numerically |16 l^Jl^ , using different meth- 
ods and reaching, to some extent, different conclusions. 

Triaxiality, anisotropy and clumping in the WIMP ve- 
locity distribution could potentially have a significant ef- 
fect on the WIMP direct detection signal. Constraints 
(and in the future possibly best fits) calculated assuming 
a standard maxwellian halo could be erroneous JTof . On 
the other hand, more optimistically, it might be possible 
to derive useful information about the local velocity dis- 
tribution, and hence the formation of the galactic halo, if 
WIMPs were detected (l6|jl|]. Belh et. al. Q have re- 
cently reanalyzed the DAMA collaborations annual mod- 
ulation signal ijl^ for a range of halo models, finding 
that the allowed region of WIMP mass-cross-section pa- 
rameter space is significantly enlarged. This illustrates 
that it is important to take into account uncertainties in 
halo modeling when comparing exclusion limits and/or 
allowed regions from different experiments. 

Given the importance of the local dark matter distri- 
bution for direct detection experiments we devote Sec. || 
to a detailed discussion of the global properties of real 
and simulated dark matter halos and recent work on the 
local dark matter distribution ]l^ , |l3| , p^ . In Sec. Ill we 
examine the effect of realistic halo modeling on exclusion 
limits. We first investigate triaxial and anisotropic halos 
models, with parameter choices motivated by the obser- 
vations and simulations, and then, more speculatively 
examine the possible effects of small subhalos. 



trary to the standard picture of halo formation in CDM cos- 
*An exception is Sikivie's late infall model Q, which con- mologies (see e.g. Refs. jl|,|T|) assumes axial symmetry and 

cold collapse. 



1 



II. GALAXY HALOS 



A. Global properties 



Observational constraints on the structure of dark mat- 
ter halos depend on the relation of luminous tracer pop- 
ulations to the underlying density distribution, and are 
complicated by galactic structure and projection effects. 
Ref. concludes that in the outskirts of spiral galaxies 
the intermediate-to-long axis ratio is likely to be greater 
than 0.8, while the short-to-long axis ratio is largely un- 
constrained with values in the range 0.3 — 1.0 reported, 
with some correlation between the method used and the 
value obtained (see Refs. and for details and 
references) . For the Milky Way (MW) , analysis of local 
stellar kinematics gives an estimate for the short-to-long 
axis value of 0.7 ±0.1 j20|, while the great circle tidal 
streams observed from the Sagittarius dwarf galaxy rule 
out ratios of less than 0.7 in the outer halo at high con- 
fidence | |2l[ | (in a flattened potential angular momentum 
is not conserved, so that orbits precess and tidal streams 
lose their coherence). The (an)isotropy of velocities in 
the MW halo is even harder to probe, however there is 
some evidence that galacticglobular clusters may have 
preferentially radial orbits p2| . 

Given the difficulties involved in 'observing' galaxy ha- 
los it makes sense to turn to numerical simulations for 
information on their possible structure. Current sim- 
ulations of galaxy halos within a cosmological context 
can resolve sub-kpc scales (see e.g. |l^j2^). Discrep- 
ancies between the halos produced in these simulations 
(which have lots of surviving dwarf galaxy sized subha- 
los and steep central profiles) and observations, have led 
to claims of a crisis for the cold dark matter model (see 
Ref. and references therein for extensive discussion). 
Most relevant for the local dark matter distribution is the 
subhalo problem which may be, at least partly, due to 
complications in comparing the observed luminous mat- 
ter with the dark matter distribution from the simula- 
tions. In particular it has been argued that gas accretion 
onto low mass halos may be inhibited after reionization 
so that a large fraction of the subhalos remain dark [ 25j . 
It has also been shown that if the observed dwarf galax- 
ies themselves have dark halos, then their masses have 
been underestimated and correcting for this would go to- 
ward resolving the discrepancy with observations p6| , ^ . 
The survival of subhalos is at least partly due to their 
concentrated profiles, so any modification to the simula- 
tions which produced halos with shallower central profiles 
could also reduce the number of surviving subhalos. De- 
spite the ongoing debate regarding the detailed compar- 
ison of the small scale properties of simulated halos with 
observations, cosmological simulations may still provide 
us with useful information about the global properties of 
dark halos. 

The shape of simulated halos varies, not just between 
different halos of the same mass, but also as a function 



of radius within a single halo, strongly if the halo has 
undergone a major merger relatively recently. Two high 
resolution Local Group halos studied in detail in Ref. |l^] 
have axis ratios of 1 : 0.78 : 0.48 and 1 : 0.45.0.38 at the 
solar radius and 1 : 0.64 : 0.40 and 1 : 0.87 : 0.67 at the 
virial radius. Adding dissipative gas to simulations tends 
to preserve the short-to-long axis ratio while increasing 
the intermediate-to-long axis ratio . 
The anisotropy parameter P(r), defined as 
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where < Vg >, < > and < > are the means of 
the squares of the velocity components evaluated at ra- 
dius r, also varies with radius. Typically f3(r) grows, 
although not monotonically, from roughly zero in the 
center of the halo to close to one at the virial radius, 
with non-negligible variation between halos (see Fig. 23 
of Ref. 1^^). The high resolution galactic mass halos 
studied in Ref. have /3{Rq) in the the range 0.1-0.4, 
corresponding to radially biased orbits. 



B. Local dark matter distribution 

In CDM cosmologies structure forms hierarchically, 
from the top down [Q. Small objects (often known as 
subhalos) form first, with larger objects being formed 
progressively via mergers and accretion. The internal 
structure of large galaxy size halos is determined by the 
dynamical processes which act on the accreted compo- 
nents. Dynamical friction causes subhalos with mass 
M > 1O^M0 to spiral toward the center of their parent 
halo within a Hubble time js^] , while the tidal field of the 
main halo can strip material away from a subhalo 26 
producing tidal streams along its orbit ||3^ . 

The local dark matter distribution can not be probed 
directly by cosmological simulations; the smallest subha- 
los resolvable in the highest resolution simulations have 
mass of order IO^A/q and it is not possible to fully re- 
solve substructure within subhalos. Little substructure 
is found within the central regions of simulated halos, 
however it is not known whether the subhalos have been 
destroyed by tidal stripping or if this is purely a resolu- 
tion effect ||l^. This is crucial for the local dark matter 
distribution as the solar radius {Rq ~ 8 kpc) is small 
compared with the radius of the MW halo, which from 
observations is thought to be in excess of 100 kpc [ p^ , 
while simulated halos with the same peak circular ve- 
locity as the MW have virial radii[] of order 200 kpc. 



^The virial radius is tlie radius which separates the virialized 
and infall regions of simulated halos and has a mean density 
within it of 178 (100) times the critical density at z = in 
the standard (A) CDM cosmology H. 
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The computing power required to directly probe the local 
dark matter distribution will probably not be available 
for a decade or so therefore other numerical [ p^Jl4| 
and semi-analytic p6[ approaches have been used to ad- 
dress the problem. 

Stiff, Widrow and Frieman |]l6| employ a semi- 
analytical approach, calculating the subhalo distribution 
as a function of mass and accretion redshift using the 
extended Press-Schechter formalism and then evolv- 
ing individual late accreting smooth subhalos within a 
growing spherical halo, to find the probability distribu- 
tion of the over-density at the solar radius. They find 
that there is a high (of order one) probability that there 
is a density enhancement of ~ 3% of the mean halo den- 
sity in the solar neighborhood, and the probability of an 
enhancement roughly equal to the background density is 
non-negligible (of order 0.01)|^. They show that if the ve- 
locity of the clump, or stream, with respect to the Earth 
is high enough it will produce a shoulder in the differen- 
tial event rate at high energies. 

Moore et. al. Q take a region with dynamical prop- 
erties similar to the Local Group, resimulated at higher 
resolution from a standard CDM cosmological simula- 
tion, and identify a subhalo with mass similar to the 
Draco dwarf galaxy. They then resimulate this subhalo 
up until its merger with the parent halo, at which point 
roughly 10% of its mass is in the form of (sub)subhalos. 
To assess the effect of tidal friction, subhalos with smooth 
density profiles are evolved within a smooth Galactic po- 
tential and it is found that small subhalos orbiting at the 
solar radius which are accreted onto the Galactic halo 
early, retain most of their mass due to their high central 
densities. 

Moore et. al. conclude that the phase space distri- 
bution at the solar radius will depend crucially on the 
Galaxy's merger history and on the internal structure of 
the smallest subhalos, arguing that it is possible that the 
local dark matter density could be zero or that a single 
dark matter stream with small velocity dispersion could 
dominate or that many tidal streams could overlap to 
give a smooth distribution. The solution to this problem 
depends on the extent to which the substructure within 
subhalos is destroyed prior to their accretion onto the 
main halo. The subhalos are of course also formed hi- 
erarchically themselves from smaller subhalos. The free- 
streaming length of ncutalinos is so small that the first 
clumps to form have mass M ~ 1O~^^M0 however, 
and to follow the accretion and destruction of subhalos 
though such a large hierarchy of scales would be decid- 
edly non-trivial. 



Helmi, White and Springel take a cosmological 
ACDM simulation, where the second largest cluster has 
been resimulated at higher resolution and then scaled 
down in size to match the MW. Most of the mass in the 
inner halo has been in place for 10 Gyr and the small- 
est resolvable subhalos do not survive to the present day 
at the solar radius. They argue that the only way a 
small subhalo with high density could exist in the solar 
neighborhood at z = is if it were first accreted by a 
large halo which subsequently had a major merger with 
the main progenitor of the Galactic halo. They conclude 
that the local dark matter velocity distribution is well 
approximated by a smooth multi-variate gaussian, with 
clumps of high velocity particles present if the MW halo 
has undergone a recent major merger. 



In summary numerical simulations produce galaxy ha- 
los which are significantly triaxial and anisotropic, with 
the shape and anisotropy of a halo depending on its in- 
dividual merger history, a picture which is broadly sup- 
ported by observations. This indicates that, even if the 
local velocity distribution is relatively smooth, the stan- 
dard spherical isotropic maxwellian halo model may not 
be a good approximation. Furthermore galaxy halos are 
formed hierarchically from the accretion of smaller sub- 
halos, which may not be completely destroyed by tidal 
friction, so that the local dark matter could be distinctly 
non-smooth. It is even possible that the dark matter 
could be distributed largely in small dense clumps. There 
is currently no consensus on the local dark matter ve- 
locity distribution however, with the results obtained 
depending on the method used to extrapolate to small 
scales below the resolution limit of cosmological simula- 
tions, so this possibility should be regarded as specula- 
tive. 



III. EFFECTS ON EXPERIMENTAL ANALYSIS 

The differential WIMP event rate due to scalar inter- 
actions can be written in terms of the WIMP scattering 
cross section on the proton, dp ]3^ : 



dR , 



Po.3 {mp + m^y 



A^T{E)F\q) 



(2) 



where the local WIMP density, is normalized to a 
fiducial value po.3 = 0.3 GeV/cm^, such that ( = Px/po.3: 
niA is the atomic mass of the target nuclei, E is the recoil 
energy of the detector nucleus, and T{E) is defined as |3^] 
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■'■These probabilities are lower limits as they only include the 
contributions of subhalos accreted after z = 1. 



where fv is the WIMP speed distribution in the rest 
frame of the detector, normalized to unity, and Wmin is 
the minimum detectable WIMP speed 
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-1.33 
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TABLE I. The values of the anisotropy parameter 7 in 
the logarithmic eUipsoidal model required to produce P — 0.1 
and 0.4. 

_ ( Eirn^ + rn^y/' (4) 



A. Triaxiality and anisotropy 

We will first examine the change in the WIMP speed 
distribution for triaxial and anisotropic, but still smooth, 
halo models. To date two self-consistent triaxial and/or 
anisotropic halo models have been studied in relation 
to WIMP direct detection: the logarithmic ellipsoidal 
model and the Osipkov-Merritt anisotropy model , 
studied in Ref. [0. We will extend the previous work 
by focusing on parameters which span the range of halo 
properties discussed in Sec. || above and including the 
detector response. 

The logarithmic ellipsoidal model ^ is the simplest 
triaxial generalization of the isothermal sphere and the 
velocity distribution can be approximated by a multi- 
variate gaussian on either the long or the intermediate 
axis (see Appendix A for further details)]^. We con- 
sider parameter values p = 0.9, q = 0.8 corresponding 
to axial ratios 1 : 0.78 : 0.48 (as used in Ref. H), and 
p = 0.72, q = 0.70 corresponding to 1 : 0.45 : 0.38, and 
locations on the long and intermediate axes. The first 
set of axial ratios is typical of the values found in sim- 
ulations and roughly consistent with observations, while 
the second is arguably rather extreme. We use values of 
the anisotropy parameter 7 (which in the spherical limit 
p = g = 1 is related to /3 by —7 = 2/3) which give /3 = 0.1 
and 0.4, and are tabulated in Tabic I. The speed distribu- 
tions, in the rest frame of the Sun normalized to unity, 
are plotted in Fig. |l| along with that for the standard 
maxwellian halo model. For both positions the triaxial 
models have a wider spread in speeds than the standard 
model, so that the differential event rate will decrease 
less rapidly with increasing recoil energy, but the change 



^Of course there is no reason to expect the Sun to be located 
on one of the axes of the halo. 




FIG. 1. The speed distributions, in the rest frame of the 
Sun, for the standard halo model (solid line), and triaxial 
models on the major axis (upper panel) and intermediate 
axis (lower panel) for p=0.9, q=0.8 and /3 = 0.1/0.4 (dot- 
ted/short dashed) and for p=0.72, q=0.7 and /? = 0.1/0.4 
(long dashed/dot dashed). 

is small on the major axis. This is because the change 
in the speed distribution is largely determined by the ve- 
locity dispersion in the direction. On the major axis, 
for parameter values which give 0.1 < f3 < 0.4, all three 
components of the velocity have roughly the same dis- 
persion, whereas on the intermediate axis the velocity 
dispersion in the </> direction is significantly larger than 
that in the z direction. Note that the speed distributions 
we consider deviate less from the standard maxwellian 
than those considered by Refs. |^ and |l^ as they use 
more extreme values for 7 (namely 16 and -1.78). 

In the Osipkov-Merritt (OM) model (see Appendix B 
for further details), which assumes a spherically sym- 
metric density profile, the velocity anisotropy varies as 
a function of radius as 

/5M = ^, (5) 

so that the degree of anisotropy increases with increas- 
ing radius, as is found in numerical simulations. Fol- 
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FIG. 2. The speed distributions for the standard 
halo model (solid line), and the OM anisotropy model 
with /3 = 0.13,0.31 and 0.4 (dotted, short-dashed, and 
long-dashed) . 



lowing Ref. ]To| we assume a NFW |9| density profileQ 
with scale radius — 20 kpc. We use values of the 
anisotropy radius ra = 20, 12, and 9.8 kpc which corre- 
spond to /3(i?o) = 0.14,0.31 and 0.4 respectively. For 
the first two values analytic fitting functions for the dis- 
tribution function have been provided by Widrow WO 
The resulting speed distributions are plotted in Fig 
along with that for the standard maxwellian halo model. 
The excess at large v is due to the increased number of 
particles on very elongated, nearly radial orbits 

To assess the effect of changes in the speed distribu- 
tion on exclusion limits we need to take into account 
the detector response, including the difference between 
the observed energy of an event and the actual recoil 
energy, non-zero energy threshold and energy resolution 
(see Ref. for further details), as these factors may 
blur out the effects of changes in the speed distribu- 
tion. We consider a Ge^^ detector with energy threshold 
Et — 4 keV with the same properties (resolution and 
form factor) as that used by the IGEX experiment 
which is optimized for detecting double-beta decay. The 
resulting differential event rates, per kg per day per keV, 
are plotted in Fig. ||, for the OM speed distributions in 
Fig. 0, for WIMPs with mass = 50 GeV and cross- 



section (Tn 



10- 



We see that, in this model, the 



differential event rate does not deviate linearly from that 
of the standard halo model as the degree of anisotropy is 
increased. Since future detectors, optimized for WIMP 
detection, will have lower thresholds and better energy 
resolution, we also plot the differential event rates for 
an 'ideal' (i.e. completely unrealistic) A = 76 detector 
where the full recoil energy is detected, the energy resolu- 
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FIG. 3. The differential event rate for the OM anisotropy 
model, with speed distributions as plotted in Fig. 2, for 



= 50 GeV and (jp 



10-*''m^for the IGEX detector (up- 
per panel) and for an 'ideal' Ge detector (lower panel). 

tion is perfect and E^^ = kcV |^. The difference between 
the differential event rates is then largest at small recoil 
energies, and would therefore be more significant for an 
experiment with a lower threshold energy. 

In Fig. ^ we plot the exclusion limits found from the 
IGEX data by requiring that the data in no more than 
one energy bin exceeds its 99.77% confidence limit, so 
as to produce 90% overall confidence limits |42j, for the 
logarithmic ellipsoidal model and for the OM anisotropy 
model. We also plot the exclusion limits from the 
Heidelberg-Moscow (HM) experiment Q| for the OM 
anisotropy model in Fig. ||. Comparing Figs. |^ and |^ 
we see that the change in the exclusion limits depends 
not only on the halo model under consideration, but also 
on the data being used; for IGEX the change in the ex- 



**Varying the inner slope of the density profile does not sig- 
nificantly affect the local velocity distribution [[10|. 



^^Somewhat counter-intuitively, at low energies the differen- 
tial event rate for the IGEX detector is higher than that for 
the ideal detector, due to the finite energy resolution of the 
IGEX detector. 
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FIG. 4. The exclusion limits from the IGEX experiment 
for the logarithmic ellipsoidal model, location on the interme- 
diate axis (upper panel) and for the OM model (lower panel) . 
Line types as in the lower panel of Fig. 1 and Fig. 2 respec- 
tively. 



elusion limits is largest for large WIMP masses, while for 
HM the change is largest for small WIMP masses. For 
different WIMP masses, different energy ranges can be 
most constraining; for the IGEX data the lowest energy 
bin is always the most constraining, while for HM as the 
WIMP mass increases the constraint comes from higher 
energy bins. It should therefore be borne in mind when 
comparing exclusion limits from different experiments, 
that changing the assumed WIMP speed distribution will 
affect the limits from different experiments differently. 

The change in the exclusion limits is not huge (of order 
tens of per-cent) for the experiments we have considered, 
however these experiments are not optimized for WIMP 
detection. As illustrated in Fig. || the change in the dif- 
ferential event rate, and hence the exclusion limit, would 
be significantly larger for an experiment with better en- 
ergy resolution and lower threshold energy. We have also 
seen that different models with the same value for the 
anisotropy parameter (3 have very different speed distri- 
butions, and hence a different effect on the exclusion lim- 
its. Furthermore it is conceivable that the local WIMP 
velocity distribution may deviate even further from the 
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FIG. 5. The exclusion limits from the HM experiment for 
the OM model. Lines as in Fig. 2. 



standard maxwellian distribution, than the models that 
we have considered. 



B. Clumps 

Even if dynamical processes produce a smooth back- 
ground dark matter distribution, late accreting clumps 
may lead to a local density enhancement and velocity 
clumping |^ 14 1, and produce a shoulder in the differen- 
tial event rate, if their density and velocity with respect 
to the earth are large enough. For the experiments we 
have been considering the lower energy bins are most 
constraining, so that only very rare high density and ve- 
locity streams would have a non-negligible effect on the 
exclusion limits. The effect of these late accreting clumps 
on the annual modulation and directional signals would 
be more significant however |p^ . 

We will now turn our attention to the consequences of 
the more speculative possibility that small subhalos may 
survive at the solar radius. We could then be located 
within a subhalo with local density in excess of the mean 
value of 0.3 GeVcm"^, on the other hand it is even pos- 
sible that we could be in a region between clumps and 

In the 



streams where the WIMP density is zero |13 



latter case all attempts at WIMP direct detection would 
be doomed to failure, and exclusion limits would tell us 
nothing about the WIMP cross-section. At the other ex- 
treme a tiny subhalo at the earth's location would pro- 
duce a distinctive signal and, due to the enhanced den- 
sity, make it easier to detect WIMPs of a given cross- 
section. Subhalos with M ^ lO^M© would have negligi- 
ble velocity dispersion and hence a delta-function speed 
distribution. The resulting theoretical differential event 
rate would be a step function with amplitude inversely 
proportional to the speed of the subhalo with respect to 
the earth, and position increasing with increasing relative 
speed and WIMP mass. Consequently for small subhalo 
velocities and WIMP masses there would be no constraint 
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on the WIMP cross-section (no WIMPs would have suf- 
ficient energy to cause an observable recoil), but as the 
WIMP mass is increased the constraints would become 
much tighter as then all the WIMPs would be energetic 
enough to cause events of a given recoil energy. In Fig. ^ 
we plot the exclusion limits on Cfip from the IGEX data 
for subhalos with various relative speed^ dominating 
the local WIMP distribution. Unlike the smooth halo 
models for this data set, as the WIMP mass is increased 
higher energy bins are most constraining, and this leads 
to sharp changes in the exclusion limits due to the sharp 
transition from no WIMPs to all WIMPs having sufficient 
to energy cause events of a given recoil energy. Note that 
the high density of a subhalo would lead to ^ > 1 , so that 
the exclusion limits on the WIMP cross-section would be 
tighter (by some unknown factor) than for the standard 
halo model. 

How dense could small scale clumps be? Dark mat- 
ter clump densities are usually parameterized in terms of 
their concentration, c, defined as 



Rv 



(6) 



where i?vir is the virial radius and rg is the scale radius, 
the radius at which the effective logarithmic slope of the 
density profile is equal to -2 (see e.g. Refs. |43p4[| ). Low 
mass halos typically form earlier, when the density of the 
Universe is higher, so that their concentrations are higher 
than those of larger halos and Bullock et. al. and 
Eke, Navarro and Steinmctz (ENS) Q have recently 
constructed toy models which reproduce the scaling of 
concentration with mass and redshift for the simulated 
galaxy-sized halos. 

Since the first neutralino clumps to form have mass of 
order lO~^^Af0 |3^, the first generation of halos which 
are formed in the same way as galaxy sized halos (via the 
accretion and merger of smaller clumps) will have mass of 
order 10~^°Mq, and virial radii of order 0.01 pc. While 
these clumps are tiny by cosmological standards, they 
are still relatively large compared to the distance trav- 
eled by the Sun in a year (10"* pc). Extrapolating the 
Bullock et. al. and ENS toy models, way beyond their 
intended range of applicability, produces significantly dif- 
ferent 'guesstimates' for c for M ^ halos of 150 
and 30 respectively [Q. Subhalos within larger halos 
may in fact be more concentrated than isolated halos of 
the same mass as dense regions tend to collapse earlier 
and tidal stripping may steepen their density profiles . 
The tidal radius |32| of a lO^^^A/© subhalo orbiting in 
the MW halo at the solar radius is of order 10"'^ pc, so 
the dense central regions of these subhalos should survive 
tidal stripping ||l^,^ . The density of simulated halos di- 
verges toward their center, so as a measure of the typical 
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FIG. 6. The exclusion limits from the IGEX data if 
a subhalo with negligible velocity distribution moving with 
speed 200, 400, 600 kms~^ relative to the earth (dotted, short 
dashed, long dashed) dominates the local density. 



density of a subhalo we use the mean density within the 
scale radius, assuming a NEW density profile ||39| ]: 



p{r) oc 



r(l -|- r /rsY 



(7) 



■'■"''In reality the subhalo speed would be another unknown 
variable. 



For c = 30(150) this gives p = 0.9(70)GeVcm-3, 
or equivalently C — 3(200). We emphasize that this 
rough calculation relies on the extrapolation of the 
concentration-mass scaling models far beyond the regime 
in which they have been numerically tested (isolated 
galaxy mass halos) and the numbers obtained should 
probably not be taken seriously. It does illustrate, how- 
ever, that i/ small subhalos survive then they may lead to 
a significant local enhancement in the WIMP density, and 
that understanding the small scale structure of galactic 
halos is crucial for realistic modeling of the WIMP veloc- 
ity distribution. 



IV. DISCUSSION 

Rapid progress is being made in the field of WIMP 
direct detection, with experiments closing in on the sen- 
sitivity required to detect neutralinos, if they constitute 
a non-negligible fraction (greater than 10"'*) of the halo 
density Q. Data analyzes usually assume the simplest 
halo model: an isothermal sphere with maxwellian ve- 
locity distribution. There is no clear justification, cither 
observational or theoretical, for this assumption apart 
from simplicity. In fact numerical simulations |p^ , p^ , ^ 
and observations ||l^,|2^,^ suggest that galaxy halos are 
triaxial and anisotropic. The local density distribution 
may also be non-smooth, with late accreting subhalos 
leading to velocity clumping |l^,|l^. More speculatively 
it is even possible that the dark matter could be mainly 
distributed in tiny dense clumps, so that the local den- 
sity distribution could be dominated by a single clump. 
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or could even be zero It is therefore crucial to ex- 
amine the effect of realistic halo modeling on the WIMP 
direct detection signal. 

In this paper we have investigated the change in ex- 
clusion limits due to triaxiality, velocity anisotropy and 
small scale clumping, taking into account detector per- 
formance and using parameter values motivated by nu- 
merical simulations and observations. Triaxiality and ve- 
locity anisotropy lead to non-negligible changes in the 
exclusion limits, even for detectors with relatively poor 
energy resolution. Furthermore the changes are different 
for different data sets and depend on how the anisotropy 
is modeled. If the local WIMP distribution is dominated 
by small scale clumps then the local density may be zero 
(making it impossible to detect WIMPs) or significantly 
enhanced (making it easier to detect WIMPs with a given 
cross-section), and the exclusion limits are changed dra- 
matically. Clearly the survival of subhalos at the solar 
radius is a very important issue for WIMP direct detec- 
tion. 

Even if the local WIMP distribution is smooth, to de- 
rive reliable constraints on WIMP parameters and com- 
pare results for different experiments a framework needs 
to be developed for dealing with the uncertainty in the 
WIMP velocity distribution; either the development of a 
framework for parameterizing deviations from a baseline 
model, or the establishment of an agreed set of bench- 
mark models, spanning the the range of plausible WIMP 
velocity distributions. 



and on the major axis 



(2 + 7)(p-2 + -7-2-l) 



2 



2(p-2 + 9-2-1) ' 
^§(2p-^"l) 

2(j5-2 + 9-2 - 1) ' 



(A3) 



where p and q are constants which satisfy < < 1 and 
are related to the axial ratios of the density distribution, 
h.2i by 



2 _ p ^jp^q^+p^ - g2) 
q2 _^ p2 _ p2q2 ' 

^2 ^ g2(p2^2 _p2 ^q2-^ 
^ (?2 ^ p2 _ p2q2 



(A4) 



and 7 is a constant isotropy parameter, which in the 
spherical limit p = q — 1 is related to f3 (as defined 
in eq. (0)) by -7 = 2/3. 



APPENDIX B: OSIPKOV-MERRITT MODEL 

The distribution function of a self-gravitating system 
with spherically symmetric density profile p{r) is given, 
assuming an isotropic velocity distribution, by Edding- 
ton's formula p2| (see also Refs. @@): 
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APPENDIX A: LOGARITHMIC ELLIPSOIDAL 
MODEL 

The logarithmic ellipsoidal model ||^ is the simplest 
triaxial generalization of the isothermal sphere and the 
velocity distribution can be approximated by a multi- 
variate gaussian on cither the long or the intermediate 
axis: 



fiv) = 



1 



(27r)3/2crr(T0Cr: 

On the intermediate axis 



exp 



(Al) 



9 —4 



(2 + 7)(l -p-2 + g-2) 

vl{2q-^-p-'^) 
2(l-p-2 + g-2) ' 

2(l-p-2 + g-2) ' 



(A2) 



1 



V87r2 



d2p d* 



1 



d*2 ^/F^ 



dp 
d* 



(Bl) 



00), <f>(r) is the potential of 
00) = -E'kin + *(r), and E 



where '^{r) = -$(r) + $(r 
the system, e = —E + ^{r 
and i?kin are the total and kinetic energy respectively. In 
the Osipkov-Merritt model |38| the distribution function 
also depends on the angular momentum, L, of the system 
through the variable Q: 



Q 



e — 



2r2 ' 



(B2) 



where ra is the anisotropy radius, which is related to (3 

by 



/3(r) 



(B3) 



so that the degree of anisotropy increases with increas- 
ing radius. The distribution function in the Osipkov- 
Merritt model is found by replacing e by Q and p{r) 
by (1 -I- r'^/rDp{r) in eq. (pl|), which then has to be 
solved numerically. Analytic fitting functions for have 
been provided by Widrow Effl for selected values of Ta, 
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for the NFW density profile |39] . Note that physical mod- 



els only exist for Ta 
potential. 



> 7'a,min! where Ta.min depends on the [22] 



[23' 
[24 
[25 



[1] L. Bergstrom, Re nt. Prog. Phys. 63 , 793 (2000). 
[2] A. Benoit et. al., |astro-ph/0206271 . 



[10 
[11 



[12: 



[is: 

[19 

[20 
[21 



[3] R. Abusaidi et. al.. Phys. Rev. Le tt. 84, 5699 (2000); D. 
Abrams et. al., 3,stro-ph/020350[ . 

[4] L. Baudis et. al., Phys. Rey D 59, 022001 (1999). 

[5] A. Morales et. al., Phys. Lett. B 532, 8 (2002). 

[6] R. Schnee, to appear in proceedings of 'COSMO 2000', 
Cheju Island, K orea, Sept. 4-8. 

[7] L. Baudis et. al., astro-ph/0005568, to appear in proceed- 
ings of "Sources and detection of dark matter in the uni- 
verse", Marina del Rey, CA, edited by D. Cline, Springer. 

[8] A. K. Drukier, K. Freese and D. N. Spergel, Phys. Rev. D 
33, 3495 (1986); F. Donato, N. Fornengo and S. Scopel, 
Astropart. Phys. 9, 247 (1998); M. Kamionkowski and A. 
Kinkhabwala, Phys. Rev. D 57, 3256 (1998); C. J. Copi, 
J. Heo and L. M. Krauss, Phys. Lett. B461, 43 (1999); 
M. Brhlik and L. Roskzkowski, Phys. Lett. B464, 303 
(1999), J. D. Vergados, Phys. Rev. Lett. 83, 3597 (1999), 
Phys. Rev. D 62, 023519 (2000); P. Belh et. al., Phys. 
Rev. D 61, 023512 (2000); A. M. Green, Phys. Rev. D 
63, 043005 (2001); C. J. Copi and L. M. Krauss, Phys. 
Rev. D 63, 04350 7 (2001); J. D. Vergados and D. Owen, 



astro- ph/0203293 



[9] N. W. Evans, C. M. Carollo and P. T. de Zeeuw, Mon. 
Not. Roy Astron. Soc. 318, 1131 (2000). 
P. UUio and M. Kamionkowski, JHEP 0103:049 (2001). 
J. D. Vergados, Phys. Rev. D 63 063511 (2001); A. M. 
Green, Phys. Rev. D 63 103003 (2001); G. Gelmini and 
P. Gondolo, Phys. Rev. D 64 023504 (2001). 
P. Sikivie, I. I. Tkachev and Y. Wang, Phys. Rev. Lett. 
75, 2911 (1995), Phys. Rev. D 56, 1863 (1997); P. Sikivie, 
Phys. Lett. B 432, 139 (1998). 
B. Moore et. al., Phys. Rev. D 64, 063508 (2001) 
A 



Helmi, S. D. M. White and V. Springel, astro- 



ph/0201289 



A. Klypin, et. al. Astrophys. J, 522, 82 (1999); B. Moore 
et. al., Mon. Not. Roy Astron. Soc. 310 1147 (1999). 
D. Stiff, L. M. Widrow and J. Frieman, Phys. Rev. D 64, 
083516 (2001). 

P. Belh, R . Cerulli, N. Forengo and S. Scopel, ^ep 



ph/0203242 



R. Bernabei et. al., Phys. Lett. B389, 757 (1996); ibid 
B408, 439 (1997); ibid B424, 195 (1998); ibid B450, 
448 (1999); ibid B480, 23 (2000). 

P. D. Sackett, Galaxy Dynamics, ASP Conf Series, 182, 
p393, eds. D. Merritt, J.A. Sellwood and M. Valluri, 
(1999). 

R. P. Oiling and M. R. Merrifield, Mon. Not. Roy. Astron. 
Soc. 311, 361 (2000). 

R. Ibata, M. Irwin, G. F. Lewis and A. Stolte, Astrophys. 
J 547, 133 (2001); R. Ibata, G. F. Lewis, M. Irwin, E. 



Totten and T. Quinn, Astrophys. J. 551, 294 (2001). 
M. Odenkirchen, P. Brosche, M. Geffert and H. J. Tu- 
cholke. New. Ast ronomy 2, 477 (19 97). 



[26: 
[27 
[28: 

[29: 

[30 
[31 

[32: 
[33 
[34: 

[35: 

[36 

[37 

[38: 

[39 

[40 
[41 

[42' 
[43^ 

[44: 

[45 

[46: 



C. Power et. a l., |astro-ph/0201544 , 
A. Tasitsiomi, istro- ph/0205464 

J. S. Bullock, A. V. Kravtsov, and D. Weinberg, Astro- 
phys. J, 539, 517 (2000). 
E. Hayashi, J. F. Navarro, J. 



Quinn, astro- ph/0203004 



E. Taylor, J. Stadel, 



F. Stoehr, S. D. M . White, G. Tormen, and V. Springel, 
astro-p h/0203342| . 



N. Katz and J. E. Gunn, Astrophys. J. 377, 365 (1991); 
J. Dubinski, Astrophys. J. 431, 617 (1994). 
T. Fukushige and M. Junichiro, Astrophys. J., 557, 533 
(2001). 

B. Moore, T. Quinn, F. Governato, J. Stadel, and G. 
Lake, Mon. Not. Roy. Astron. Soc. 310, 1147 (1999). 
S. D. M. White and M. Rees, Mon. Not. Roy. Astron. 
Soc. 183, 341 (1978); G. R. Blumenthal, S. M. Faber, J. 
R. Primack and M. J. Rees, Nature 311, 517 (1984); M. 
Davis, G. Efstathiou, C. S. Frenk and S. D. M. White, 
Astrophys. J. 292, 517 (1985). 

J. Birmey and S. Tremaine, Galactic Dynamics (Prince- 
ton University Press, Princeton, NJ 1987). 
A. Helmi and S. D. M. White, Mon. Not. Roy. Astron. 
Soc. 307, 495 (1999). 

S. Cole and C. Lacey, Mon. Not. Roy. Astron. Soc. 281, 
716 (1996); V. R. Eke, S. Cole and C. S. Frenk, Mon. 
Not. Roy Astron. Soc. 282, 263 (1996); G. Bryan and 
M. Norman, Astrophys. J. 321, 80 (1998). 
J. R. Bond, S. Cole, G. Efstathiou, and N. Kaiser, As- 
trophys. J. 379, 440 (1991); C. Lacey and S. Cole, Mon. 
Not. Roy Astron. Soc. 262, 627 (1993). 
S. Hofmann, D. J. Schwarz and H. Stocker, Phys. Rev. 
D 64, 083507 (2001). 

G. Jungman, M. Kamionkowski and K. Griest, Phys. 
Rep. 267, 195 (1996). 

L. P. Osipkov, Pis'ma Astron. 55, 77 (1979); D. Merritt, 
Astron. J. 90, 1027, (1985). 

J. F. Navarro, C. S. Frenk, and S. D. M. White, Astro- 
phys. J, 462, 563 (1996); ibid 490, 493 (1997). 
L. M. Widrow, Astrophys. J. Suppl. 131, 39 (2000). 
J. D. Lewin and P. F. Smith, Astropart. Phys. 6, 87 
(1996). 

A. M. Green, Phys. Rev. D 65, 023520, (2002). 

J. S. Bullock et. al. Mon. Not. Roy. Astron. Soc. 321, 

559 (2001). 

V. R. Eke, J. F. Navarro, and M. Steinmetz, Astrophys. 
J. 554, 114 (2001). 

P. UUio, L. Bergstrom, J. Edsjo and C. Lacey, astro- 



ph/0207125 



G. Duda, G. Gelmini and P. Gondolo, Phys. Lett. B 529, 
187 (2002). 



9 



